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Tipping  Point  of  NG 

A  minority  of  committed  agents  can  persuade 
the  whole  network  to  a  global  consensus. 


The  critical  value  for  phase  transition  is  called 


Saddle  node  bifurcation 


Below  Critical 


Above  Critical 


Meanfield  Assumption  and  Complete 
Network 


The  network  structure  is  ignored.  Every 
node  is  only  affected  by  the  meanfield. 

The  meanfield  depends  only  on  the 
fractions(or  numbers)  of  all  types  of  nodes. 

Describe  the  dynamics  by  an  equation  of 
the  meanfield  (macrostate). 


Scale  of  consensus  time  on  complet 
graph 


2  word  Naming  Game  on  complete  graph 


Expected  Time  Spend  on  Each 
Macrostate  before  Consensus  (without 

committed  agents) 


n, 


n 


B 


0  0 


NG  with  Committed  Agents 

<\  q=0.06<qc  *’  q=0.12>qc 


q  is  the  fraction  of  agents  committed  in  A. 

When  q  is  below  a  critical  value  qc,  the  process 
may  stuck  in  a  meta-stable  state  for  a  very 
long  time. 


2  Word  Naming  Game  as  a 
2D  random  walk 


P(A-)  =  nA{N  -nA  +  nB)/2N2 
n_B  P(A+)  =  (N  -  nA  -  nB)(N  +  nA  -  nB)/2N2 

P(B-)  =  nB(N  +  nA-  nB)/2N2 

P(B+ )  =  (N  —  nA  —  nB)(N  -  nA  +  nB)/2N2 

Po  —  (nA  +  nB)/2N  +  ( nA  -  nB)2 /2N2 


P(B+) 


0 

Transient  State 

-> 

• 

Absorbing  State 

n  A 

Linear  Solver  for  2-Name  NG 


r(n  4,  ns)  is  the  absorbing  time  starting  from  network  state  (tia,  Ub)  . 
t(nA:  ub)  is  the  average  time  stay  in  network  state  (ua^b)  before  leaving. 

Have  equations: 

P(A+)t {nA  +  1,  ub)  +  P(A-)t(tia  -  Ijib)  +  P{B+)r(nA,nB  +  1)  +  P(B-)T[nA,nB  -  1) 
T(nA,nB)  =  - - - - - - — - - — - - - - — — - - 

+  t(nAlnB) 

Then  we  assign  an  order  to  the  coordinates,  make  r(nA,  nB)  ,  t , ns ) 
into  vectors,  and  finally  write  equations  in  the  linear  system  form: 


f  =  Mr  +  t 


SDE  models  for  NG, 


dX.  =  fidt  + 


NG  and  NG 


Diffusion  vs  Drift 


Diffusion  scales  are  clear  from  broadening 
of  trajectories  bundles 

Drift  governed  by  mean  field  nonlinear 
ODEs  can  be  seen  from  the  average  / 
midlines  of  bundles 


Assume  a  very  natural  social  -  political  condition,  generalizing  NG 

Where  the  network  is  divided  into  k+1  sub-populations 

Each  with  a  different  fixed  propensity  to  signal/vote/utter  the  opinion  A 

When  hearing  the  word  A  a  node  from  subgroup  j  <  k  will  move  to  subgroup  j+1 
Likewise  hearing  B  a  node  from  j  >  0  will  move  to  subgroup  j-1 

The  probability  that  a  node  from  subgroup  j  will  signal  A  is  j/k 

Same  node  has  probability  1  -  j/k  of  signalling  word  B 


Additional  rules  of  k-NG 


A  node  s  is  drawn  at  random  from  the 
network  and  sends  out  a  signal  A  with  p 

or  s  signals  B  with  prob  =  1  -  p 

Next  a  node  L  is  drawn  at  random  to 
receive  the  signal  A 

or  node  L  is  drawn  to  receive  the  signal  B 


k-NG 


The  subgroups  j  =  0,  k  correspond  to 
those  nodes  that  are  completely  convinced 
of  the  B  and  A  opinion  resp.  OR 

Equivalently  those  nodes  that  with  prob  1 
Signals  B,  A  resp. 


Voting  or  Polling 


The  key  network  quantity  is  an  average 
network  opinion  obtained  by  polling: 

p  =  sum  over  subgrps  j  =  0  to  k 
of  n(j)j/kN 

It  gives  the  probability  of  a  speaker  chosen 
at  random  signalling  the  word  A 


Stochastic  Dynamics 

Derivation  of  k-dim  coupled  Random  Walk: 


For  j  =  1  to  k-1, 

n(j,t+1)  =  n(j,t)  +  1,-1  with  resp.  prob. 
P(+1)  =  p(t)  n(j-1,t)/N  +  (1-p(t))  n(j+1  ,t)/N 
P(-1)=  (1-p(t))  n(j,t)/N  +  p(t)  n(j,t)/N 
=  n(j,t)/N 


Random  walk 


For  the  distinguished  subgrps  at  both  ends 
of  the  opinion  spectrum,  j  =  0  ,  k 

n(0,t+1)  =  n(0,t)  +  1, 0,  -1  with  prob. 
P(+1)  =  (1-p(t))  n(1,t)/N 
P(0)  =  (1-p(t))  n(0,t)/N 
P(-1)=  p(t)n(0,t)/N 


Subgroup  k 

n(k,t+1 )  =  n(k,t)  +1,0, 

P(+1)  =  p(t)  n(k-1  ,t)/N 
P(0)  =  p(t)  n(k,t)/N 
P(-1)  =  (1-p(t))n(k,t)/N 


-1  with  prob. 


Shadow  Walk 


Define  a  scalar  nonlazy  random  walk  on 
the  quantity  p  whose  consensus  times  are 
lower  bounds  for  the  actual  expected 
times. 

p(t+1)  =  p(t)  +  1/kN,  -1/kN  with  prob 


P(+1)  =  p(t) 

P(-1)=  1-P(t) 


SDE  or  Diffusion  Model 


In  continuous  rime  the  SDK  is 


d!)  =  M{l\)dt  +  Y(l\)d\\ 
M  =  (2M  -  l)/fe 


Solution  of  SDE 


now  done  using  standard  Kolmogorov  Backwards  Equation  method 


V’ 

0,  u(l,  t)  =  1  whervu(p,  t)  =  Pr{  !\  = 


«(/')  =  lim  «(/>./) 

t—'O 0 


Expected  Times  to  Consensus 


First  calculate 


=  exp  ( 2  v  N 


Next  we  get 


|0P  G{x)<k  _  Jo  exl>  (2 y/Ny/x-x*)  dx 
Jo  G{x)dx  ex])  (I'/N'fx  -  i2j  dx 


Expected  times  -  exit  times  -  stop  times 


and  the  expected  time  to  consensus  of  the  A  opinion  without  committed 
agents  and  conditioned  on  eventual  fixation  of  the  A  opinion  is  given  by 


where 


is  the  unconditioned  expected  time  to  consensus  given  that  Po  =  p  and  solves 
the  stationary  form  of  the  KBE 


<fT  2M  dT  2  u(p) 
+1 r~dp  +  ^T 

with  boundary  conditions 


liin 

p— 0 


t  <  oc 

0. 


continued 


The  solution  is  given  by 


■ 


y-y* 


Higher  stubbornness  -  same  qualitative, 
robust  result 


2 


40Q  SOD 


400 

35a 

300 

25a 

200 

150 

100 

so 

0 

0 


522 

T 


3DD 


100 


<  =  3,  N  = 


200 


urtrtomi,  p=Q.S 
centered.  p=0.5 
5'  zed  p=2  5 
centered,  p-0.6 
jc  anzed.  p=0.4 


Figure  l :  Simulations  of  the  Naming  Game  with  K=3  and  N  =  500,  Note  the 
existence  of  a  center  manifold  to  which  all  trajectories  tend.  The  simulations 
that  start  far  away  from  the  center  manifold  approach  it  lief  ore  drifting  to  a 
consensus  state. 


Higher  stubbornness  -  same  qualitative, 
robust  result 
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Figure  l :  Simulations  of  the  Naming  Game  with  K=3  and  N  =  500,  Note  the 
existence  of  a  center  manifold  to  which  all  trajectories  tend.  The  simulations 
that  start  far  away  from  the  center  manifold  approach  it  lief  ore  drifting  to  a 
consensus  state. 


Other  NG  variants  -  same  1 D  manifold 


No 


Figure  2:  Simulations  of  the  Naming  Game  with  K=10  and  N  =  50Q.  Some  tra- 
jpctoricH  ■  n<  ►!  Mpproarh  tin-  nut n  manifold  aw  rjuirkly  aw  they  did  in  the  1\  ‘t 
case.  Those  simulations  that  start  with  p=0. 5  approach  the  center  manifold  he- 
fore  drifting  to  the  consensus  state.  However,  the  simulations  that  start  with  a 
lews  centralized  p  value  drift  to  the  cordons uw  state  more1  quickly.  Note  that  the 
simulations  starting  in  the  polarized  state  with  p  =  U.4  come  close  to  the  center 
manifold  he  ion*  drilling  tn  confimsUS-  On  the  other  Imiid.  the  ^iim  dnt  ions  in  the 
centralized  state  with  p=0.<>  drift  to  consensus  before  any  significant  number  of 
agents  enter  .Vq, 


3D  plot  of  trajectory  bundles  - 
stubbornness  K  =  10  as  example  of 

_ _ 1 _ i  /XX  T _ ■! _ _  _ _l  XN  I  1 _  XX  XX  Jk  xx  \ 


K  =  IQ,  N  =  500 


uniform,  p;  =  0  5 
centered,  Pj  =  0.5 
pa  anzed.  p:  =  0  5 
centered,  =  0.6 
pa  anzed.  p|  =  0  4 


Figure  l:  A  3D  plot  of  the  simulations  of  the  Naming  Game  with  K=3  and  N 
=  500,  Here,  all  simulations  approach  the  center  manifold  before  p  deviates  too 
much  from  its  initial  value. 


Consensus  time  distribution 


Recursive  relationship  of  P(X,  T),  the  probability 
for  consensus  at  T  starting  from  X,  Q  is  the 
transition  matrix. 


P{nA,nB,T  +  1)  =  13(714  +  lfnB\nA,nB)P{nA  +  l,nBtT)  +  <3(^4  -  l,nB\nA}nB)P(nA  -  l,nBiT) 

+Q{nA,nB  +  l\nA7nB)P(nA,nB  -|-  1}T)  +  Q{nA,nB  -  l\nA,nB)P(nA,nB  -  1,T) 
+Q{nA>nB  +  2\nA, nB)P{nA, nB  +  2.T)  -\-Q{ua  +  %nB\nAinB)P(nA  +  %nB,T) 

Take  each  column  for  the  same  T  as  a  vector: 

P(T+l)  =  Q*P(T) 

Calculate  the  whole  table  P(X,T)  iteratively. 

Take  each  row  for  the  same  X  as  a  vector: 

nA,nB)(Tc  =  T )  =  {P(nA,nB:.T)). 


Consensus  time  distribution 


Tc  X  104  TC  X  107 


Red  lines  are  calculated  through  the  recursive  equation. 

Blue  lines  are  statistics  of  consensus  times  from  numerical  simulation(very  expensive), 
(done  by  Jerry  Xie) 


P(T)  *  std(T) 


p=0.12 


Consensus  Time  distribution 


p=0.08 


(  T-E[Tc]  )  /  std(Tc) 

Below  critical,  consensus  time  distribution  tends  to  exponential. 
Above  critical,  consensus  time  distribution  tends  to  Gaussian. 

For  large  enough  system,  only  the  mean  and  the  variance  of  the 
consensus  time  is  needed. 


NG  on  RGG 
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travelling  wave  solution  of  u(x)  when  q=0.05<qc 
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NG  on  RGG  past  Tipping  point 
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Homogeneous  Pairwis  Assumption 

The  mean  field  is  not  uniform  but  varies  for  the  nodes 
with  different  opinion. 


Make  it  rigorous: 

suppose  three  nodes  are  linked  as  i-2-3 
Xj  is  the  opinion  of  node  i 
ki  is  the  degree  of  node  i 

PfalX*)  ?  PiXJ 
P{X1\X2;XZ)  =  P{X1\X2) 

Elh  |Xj  =<k> 

P(X1\k1)  =  P(X1) 

P(X1\X1:kl,k2)  =  P(X1\X2) 


Numerical  comparison 


Trajectories  mapped  to  2D  macrostate 


0.95 


Concentration  of  the  consensus  time 


Change  of  the  tipping  point  w.r.t.  the 
average  degree 


(nB 


N=10000 


*  MeanField 
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The  local  mean  field  for  the  node  with  opinion  C : 
P(-|C)  =  [P{A\C),P{D\C),P{AB\C)\r  ,C  —  A,B,AB 
The  number  of  different  type  of  links: 

L  =  [La-a?La-b-.La-abiLb-b,Lb_ab?LAb-a  bt 


_ ,  /  P{A\A)  \ 

W)(L)  =  P(B\A) 
\P(AB\A)J 


1 

'2La-a  +  La-  b  +  La-ab 


2£j4-j4\ 
La-b  I 
La-ab } 


P(-\B)(L) 


(  P(AB)\ 
P(BB) 
\P(AB\B)j 


1 

La-b  +  2Lb~b  H-  Lq-ab 


Ba-b  \ 

2Lb_b  I 

\Lb-ab } 


P(:\AB)(L) 


(  P(A\AB)  \ 
P(B\AB) 
\P(AB\AB)J 


I 

La-ab  +  Lb-ab  +  2Lab~ab 


(Ba-ab 
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2-Bab-ab 


Analyze  the  dynamics 


1  .Choosing  one  type  of  links,  say  A-B,  and  A  is  the  listener. 

2.  Direct  change:  A-B  changes  into  AB-B. 

3.  Related  changes:  since  A  changes  into  AB,  <k>-1  related 
links  C-A  change  into  C-AB.  The  probability  distribution  of  C 
is  the  local  mean  field  P(  |A). 


Direct 


Local  mean  field  equation 


E[AL\L]  =  jj  [D  +  (<  k  >  -1)R]  L 
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R  =  ( +  QBP(.|S)],Q4[|P('M)  -  7?(-  |AS)],0,Qb[jP(-|5)  -  -aP(-\AB)\,-(Qa  +  Qb)P(-\AB)  ) 
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Normalized  equation: 


1 


<k> 


D  +  { 


<  k  >  —  1 
<k> 


<k> 


l 


SDE  model  of  NG 


r+AT-i 

AX  (AT)  =  J2  Xn(X(i)).  ^ 


i=T 


E[An] 


E(A  n) 


Ndt  dX  =  ndt  +  —r LdW 
~~  y/N 


[dX.dX] 


N 


Merits  of  SDE  model 


dX  =  j. Xdt  +  ~^=dW 

y/N 

Include  all  types  of  NG  and  other 
communication  models  in  one  framework 
and  distinguish  them  by  two  parameters. 

Present  the  effect  of  system  size  explicitly. 

Collapse  complicated  dynamics  into  1-d 
SDE  equation  on  the  center  manifold. 


Thanks 


